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The  degradation  due  to  random  segment  jitter  and  phase  mismatch 
on  the  average  far  field,  on  axis  irradiance  of  a  laser  beam  trans¬ 
mitted  from  a  telescope  with  a  segmented  primary  mirror  is  computed. 
The  derived  expression  is  exact  within  the  stated  assumptions  and 
offers  a  more  accurate  alternative,  for  this  application,  to  the 
commonly  used  Strehl  approximation, 
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For  3  system  designed  to  transmit  a  laser  beam  to  a  receiver  or  target, 
there  are  four  options  available  to  increase  the  intensity  of  the  deposited 
radiation: 

1.  Use  a  shorter  wavelength  laser. 

2.  Minimize  system  errors. 

3.  Increase  laser  power. 

4.  Increase  the  size  of  the  transmitting  optics. 

Each  evolutionary  step  taken  in  advancing  the  technology  usually  incorporates 
improvements  in  at  least  the  last  three  areas.  Although  building  a  bigger  trans¬ 
mitting  telescope  might  appear  to  be  one  of  the  most  straightforward  approaches 
to  improving  system  performance ,  for  applications  involving  very  long  propagation 
distances,  it  will  become  increasingly  difficult  to  fabricate  and  control  optics 
of  the  required  size. 

An  approach  which  would  potentially  alleviate  some  of  the  problems  inherent 
in  simply  scaling  up  a  transmitting  telescope  of  conventional  design  is  to  divide 
the  laser  beam  up  among  several  smaller  telescopes  or  to  make  the  primary  mirror 
segmented.  The  advantages  of  this  approach  to  telescope  design  for  applications 
in  high  resolution  imagery  is  widely  appreciated.  The  Multiple  Mirror  Telescope 
currently  in  operation  on  Mt  Hopkins  collects  light  with  six  hexagonal  primary 
mirrors.  Proposed  plans  for  a  larger  follow-on  to  the  Space  Telescope  include 
a  segmented  primary. 
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In  either  case,  the  individual  telescopes  or  segments  would  have  to  be 
phased  together  so  that  the  pieces  of  the  laser  beam  would  add  coherently  on 
the  target.  In  addition  to  the  phasing  requirement,  the  different  subapertures 
of  the  transmitter  would  all  have  to  be  pointed  to  the  same  spot  on  the  target. 

In  practice,  these  requirements  would  be  satisfied  only  imperfectly.  There 
would  be  some  residual  phasing  error  and  some  pointing  jitter  of  the  mirror  seg¬ 
ments  which  would  vary  from  one  segment  to  the  next.  The  effect  of  these  residual 
errors  on  the  far  field,  on  axis  irradiance  is  computed  below. 

Consider  a  transmitting  aperture  of  area  A  divided  into  N  circular  subapertures 
of  area  A/N  =  nD^/A.  If  monochromatic  radiation  of  wavelength  A  is  transmitted 
through  the  N  subapertures  and  brought  to  a  focus  a  distance  z  from  the  telescope, 
then  the  field  in  the  focal  plane  is  given  by  (neglecting  purely  phase  terras)^ 


U(7)  =  (l/Az)  /  U(r.)  exp(-ikr*r./z)  d^r^ 
A 


N  _ 

=  (1/Xz)  X)  S  U  (r.)  exp(-ikr*r  /z)  d^r 

n=l  A  “ 
n 


(1) 


where  k  =  2n/X,  and  is  the  field  in  the  n—  subaperture  A^.  We  assume  that 

each  U  is  a  uniform  amplitude  wave  with  a  phase  consisting  of  random  piston 
n 

and  tilt. 


U^(r^)  =  (P/A)^  exp  |ik  [a^  +  b^  (x^ 


c 

n 


(^1 


(2) 
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where  (x  ,y  )  is  the  centroid  of  the  n—  subaperture,  and  P  is  the  total  power 


n"n 


in  the  beam.  The  irradiance  pattern  in  the  focal  plane  is  then  given  by 


I(r)  =  U(r)  U^Cr) 


=  (PA2z2A) 


n,m  A 


m 


X  exp  {ik[a^  +  b^Cxj-x^)  + 

X  exp  {*ik[a^  +  b^Cx^-x^)  +  c^(y2-yB,)H 


(3) 


We  assume  that  a^,  b^,  and  c^  are  pairwise  independent  Gaussian  random 
variables  with  mean  0  and  variances  given  by 


<  a^  > 

n 

<  b^  > 

n 


=  rr2 


<  c^  > 
n 


=  ol 


(4) 


<->  denotes  an  ensemble  average.  The  average  focal  plane  irradiance  distribution 


is 


<  I(r)  >  =  (P/\2z2a) 


I  /J  d^rj  d^r^  exp[-ikr •(r^-r2)/z] 


n  A 


X  <  exp[ikb^(Xj-X2)]  >  <  exp[ikc^(y^-y2)]  > 
+  J  d^r^  /  d^r^  exp[-ikT* (rj-r^)/z] 


n#m  A 


m 


X  <  exp(ika^)  >  <  exp(-ika^)  >  <  exp[ikb^(Xj-x^)l  > 

X  <  exp[-ikb^(x2-Xjjj)l  >  <  exp[ikc^(yj-y^)]  >  <  exp[-ikCj^(y2-y^)l  > 
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(PA2z2A) 


1  JS  exp["ikr*(rj-r2)/zl 

X  exp  {-k^  a|  [(x^-x^)^  +  (yj"y2)^l/2} 

+  ^  X  exp[-ikF*(r^-r^)/2]  (5) 

n^m  A  A 

n  m 

X  expC-k^ap  exp  {-k202 [ (Xj-x^)^  +  (y^-y^)"  + 

Transforming  to  sum  and  difference  coordinates  in  the  first  set  of  double 
integrals  and  translating  the  origin  to  the  suhaperture  centroids  in  the 
second  set  of  integrals  in  equation  (5),  it  is  straightforward  to  show  that 


<  I(r)  >  =  (P/X2z2A) 


4AD2  X  pdp  exp(-k2a2D2p2/2) 
0 


X  [cos'  (p)  -pVl-p^l  J^(kDrp/z) 


+  An^  exp(-k2ap  exp[-ikr ‘(r^-r^)] 


n^tra 


D/2 

X  [  X  Tjdr^  exp(-k2o2r2/2)  J^Ckrr^/z)]' 
Evaluating  <  IC?)  ^  on  axis,  we  obtain 


<  1(0)  >  =  (P/X2z2a) 


4AD2  X  pdp 
0 


X  exp(-k2a|D2p2/2)  [cos  (p)  -pVl“P^l 

D/2 

+  4n2  exp(-k2a2)(N2-N)  [  f  r^  exp(-k2a2r2/2)drj]= 


(6) 


(7) 
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The  integrations  appearing  in  equation  (7)  can  be  performed  analytically.  The 
first  integral  can  be  evaluated  by  integrating  by  parts  and  then  using  a  trig" 
onoBsetric  substitution  to  transform  it  to  a  standard  integral  identity  for  Bessel 
functions.  The  result  is 


<  1(0)  >  =  (PAA^z^) 


[8/(Nk2a2D2)] 


x[i  -  exp(-k2a2D2/4)  [Ij,(k2a2D2/A)  +  Ij(k2a2D2/4)]] 


+  (64/(N2k^ajD‘*)]  (N2-N)  exp(-k2aj)  (8) 

X  [1  -  exp(-k2o2D2/g)]2 

» 

Iq  and  denote  the  modified  Bessel  functions  of  order  0  and  1  respectively. 2 

Let  p^  =  a^/K  and  p  =  a^/K.  Then  p^  is  the  rms  piston  error  in  waves  and 
p  is  the  rms  tilt  error  measured  in  units  of  X/D,  the  diffraction  angle  of  the 
subapertures.  If  we  normalize  <  1(0)  >  by  its  diffraction  limited  value,  , 

we  obtain  the  Strehl  ratio  which  we  denote  by  ®  function  of  p^  and  p, 

is  given  by 

=  (2/Nn2p2)  {1  -  exp(-n2p2)[i^(n2p2)  +  Ij^(n2p2)l} 

.+  4  [(N-1)/(N7i‘*p^)]  exp(-4n2pp  [2  -  exp(-n2p2/2)12  (9) 
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Equation  (9)  can  be  compared  to  the  Strehl  approximation: 

Irei  =  exp[-4n2(N-l)p^2/N  -n^p^/2]  (10) 

The  Strehl  approximation  gives  lower  values  than  equation  (9),  but  it  becomes 
more  accurate  as  the  number  of  subapertures  increases.  For  N  S  6,  the  Strehl 
approximation  is  within  a  few  percent  of  the  correct  value  for  ^0.4. 

It  is  instructive  to  note  the  behavior  of  ior  certain  special  cases. 

Case  1.  N  =  1.  As  expected,  is  independent  of  when  there  is 

only  one  aperture. 

=  (2/7t2p2)  j  1  -  exp(-n2p2)[lQ(n2p2)  +  Ij(7t2p2)]  |  (11) 

Case  2.  =  0.  The  limit  as  p  ->  0  can  be  evaluated  by  substituting  the 

power  series  for  the  exponential  function  and  the  Bessel  functions  in  equation  (9) . 

Irei  =  n  +  (N-1)  exp(-4n2p^2))/N  (12) 

for  p  =  0.  Note  that  I  ,->lasp  ^0,  and  I  ,  1/N  as  p 

rel  *^o  ’  rel  *^o 

Case  3.  N  ^  For  large  N, 

=  4  exp(-4n2p^2)  Jj.  eXp(-7l2p2/2)]2/(7t4p4)  (13) 
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Care  should  be  exercised  in  interpreting  equation  (13).  If  A  is  held  constant 
as  N  then  D  0.  Then  for  a  fixed  value  of  p,  a.j,  •»  ».  In  applying  equa¬ 

tion  (13),  it  is  necessary  to  keep  in  mind  what  is  varying  and  what  is  being 
held  constant. 

If  we  let  p  0  in  equation  (13),  then 

I  ,  ->■  exp(-4n2p^2)  (14) 

Equation  (14)  is  simply  the  usual  Strehl  approximation  which  is  exact  in  the 
limit  of  phase  aberrations  with  zero  correlation  length.  That  is  just  the 
situation  represented  by  uncorrelated  piston  errors  over  N  subapertures  when 

N  ■>  ». 

V 

Figures  1-5  show  the  dependence  of  1^.^]^  P  different  values 

of  N.  In  Figure  1,  I^^j^  is  shown  as  a  function  of  p  for  the  case  N  =  1.  There 
is  no  dependence  on  p^,  and  the  simplified  expression  given  in  equation  (11)  holds. 

Figures  2,  3,  and  4,  show  results  for  the  cases  N  =  6,  18,  and  36,  respectively. 

(These  values  correspond  to  hexagonal  arrays  of  circular  apertures  with  a  hole  in 

the  center.)  In  each  figure,  1^.^^^  is  plotted  as  a  function  of  p.  The  different 

curves  in  each  figure  correspond  to  different  values  of  p^.  p^  is  varied  from  0.0 

to  0.5  in  increments  of  0.05;  i.e.,  the  piston  error  is  incremented  by  1/20  wave. 

Of  course.  I  ,  decreases  as  p  increases  so  the  top  curve  in  each  figure  corres- 
rel  o 

ponds  to  p  =0,  and  each  successively  lower  curve  corresponds  to  the  next  higher 

value  of  p  . 

o 
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Figure  5  corresponds  to  the  limit  N  ■♦  ».  The  carves  were  drawn  using  equa¬ 
tion  (13).  The  piston  error  is  again  varied  in  increments  of  1/20  wave. 

Although  many  important  sources  of  aberrations  have  been  ignored  and  some 
simplifying  assumptions  have  been  used,  equation  (9)  should  be  useful  for 
determining  the  phasing  and  pointing  accuracy  requirements  for  the  component 
pieces  of  a  multiple  mirror  or  segmented  mirror  telescope. 
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